Experimental test of quantum non-Gaussianity of heralded single photon state 
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We report on experimental verification of quantum non-Gaussianity of a lieralded single photon 
state with positive Wigner function. We unambiguously demonstrate that the generated state cannot 
be expressed as a mixture of Gaussian states. A sufficient information to witness the quantum non- 
Gaussianity is obtained from a standard photon anti-correlation measurement. 
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Quantum properties of light are exemplified by sta- 
tistical behaviors which do not admit explanation based 
on a semiclassical theory. Since coherent states repre- 
sent quantum analogue of classical coherent light, a state 
that cannot be expressed as a convex mixture of coher- 
ent states is commonly considered to be nonclassical [l| . 
In particular, during recent decades nonclassical squeezed 
states of light have become a crucial resource for quantum 
optics, metrology, and quantum information processing 

m- 

Pure squeezed coherent states represent extremal 
points of a convex set of stochastic mixtures of Gaussian 
states. All such states possess positive Wigner function 
and can be obtained from coherent laser beams using 
classical mixing and quantum interactions described by 
quadratic Hamiltonians. Exploiting higher-order nonlin- 
earities involved in the photon detection process, states 
with negative Wigner function can be conditionally gen- 
erated from the squeezed states y-[8| . Wigner function of 
these highly non-classical states exhibits a distinctly non- 
Gaussian shape that cannot be obtained as a stochastic 
mixture of Gaussian functions. 

The famous Hudson theorem establishes an equiv- 
alence between the non-Gaussianity and negativity of 
Wigner function for pure states [9|. However, this re- 
lation does not simply extend to mixed states [10[. Pre- 
vious approaches towards the non-Gaussianity witness 
or measure for mixed states do not distinguish non- 
Gaussianity which is compatible with a simple mixture 
of Gaussian states, and they also require complete infor- 



mation about the quantum state llHl3|. This brings a 
very basic and fundamental physical problem to our at- 
tention: Which mixed non-classical quantum states with 
positive non-Gaussian Wigner function do not admit ex- 
planation based solely on stochastic non-Gaussianity? 
Mathematically, we search for a witness certifying that 
the state cannot be constructed as a mixture of Gaus- 
sian states. Very recently, a directly measurable witness 
of the quantum non-Gaussianity has been theoretically 
proposed IJ]. The witness is based on knowledge of 



probabilities of vacuum and single-photon states only, yet 
it can detect a wide class of states with positive Wigner 
function which are not mixtures of Gaussian states. 

A heralded single photon source is an excellent ex- 



ample for testing the power of this witness in a labo- 
ratory. In the absence of background noise, the gener- 
ated state would be a mixture of a single photon state 
and a vacuum due to losses and imperfect coupling and 
mode- matching. If the probability of vacuum dominates, 
then the state exhibits positive Wigner function. Nev- 
ertheless, the witness [ij] still proves that it is not a 
mixture of Gaussian states. Here we apply the witness 
to approximate single-photon states conditionally gener- 
ated by detection of an idler photon from a photon pair 
produced by the process of spontaneous parametric fre- 
quency down-conversion (PDC). Our detection scheme 
consisting of a beam splitter and two single-photon de- 
tectors is the one commonly employed to test the anti- 
correlation properties of single-photon sources [l5| . This 
measurement allows us to obtain suitable estimates of 
vacuum and single-photon probabilities, which are re- 
quired for the non-Gaussianity witness. The verification 
of quantum non-Gaussianity thus conveniently comple- 
ments other typically performed non-classicality tests of 
single photon sources. 

Theory: Let Q denote the set of all mixtures of Gaus- 
sian states. We would like to show that a given state 
p cannot be expressed as a convex mixture of Gaussian 
states, p ^ Q, even though p possesses a positive Wigner 
function. This can be accomplished using a criterion re- 
cently derived in Ref. [l^. This criterion is based on 
photon number probabilities and can be expressed as an 
upper bound on single-photon probability pi for a given 
vacuum state probability p^. If the measured pi exceeds 
this bound, then p ^ Q- The bound can be derived by 
maximizing pi for a fixed po over all pure Gaussian states 
[14| and can be conveniently expressed in a parametric 
form, 

g-d^[l-tanh(r)] ^2 g-d^[l-tanh(r)] 

cosh(r) cosh'^(r) 

Here r > is the squeezing constant and the displace- 
ment reads tP = (e^"" — 1)/4. All probability pairs (poiPi) 
achievable by mixtures of Gaussian states form a convex 
set that is shown in Fig. 1(a) as a blue area. Note that 
the boundary of this area is specified by the formula ([T|) . 
In analogy with entanglement witnesses [16|, we can 



define a non-Gaussianity witness |14| . 
W{a) =ap„+pi. 



(2) 



If W{a) > Wcia) then p (f G. Tfie bound Waia) = 
maxpgg W{a) can be obtained by solving tlie equation 
(1 + e^'')a = e^'"(3 — e^'') with respect to r, determining 
Po and pi from Eq. ([l} and inserting them into Eq. ^ . 
As indicated by a green dashed line in Fig. 1(a), each line 
apo +pi = Wcia) is a tangent to the boundary curve dJ) 
and divides the plane into two half-planes. All points 
iPoiPi) lying in the half-plane apo + pi > Wcia) are 
certified by the witness to correspond to a state p ^ Q- 

Let us now consider practical determination of the 
probabilities po and pi. Since the currently commonly 
available avalanche photodiodes (APDs) are not capable 
of resolving the number of photons, one needs to em- 
ploy an advanced photon-number resolving detector |17l - 
[l9[ with demanding operation conditions or some sort 

20|427|. Perhaps the conceptu- 
shown in Fig. 1(b), is based on 



of multiplexed detector 
ally simplest scheme. 



splitting the incoming signal on a balanced beam splitter 
BS and placing an APD on each output port of the BS. 
This setup is commonly used for measurement of the g'-^^ 

factor [illiili^. 

Although this scheme provides some information about 
photon statistics there are several factors that need to be 
carefully considered. One important issue is the detec- 
tor efficiency and other losses which combine to overall 
efficiency rj. Compensation of rj would require its precise 
calibration, which is a non-trivial task. However, we can 
simply include losses into state preparation. Let i^ de- 
note a lossy channel with transmittance rj. This channel 
maps Gaussian states onto Gaussian states. Therefore, 
ii p E G then also p,, = L,^[p] g G- This implies that if 
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FIG. 1: (Color online) (a) Inferring non-Gaussianity from 
photon number probabilities po and pi of vacuum and single- 
photon states. The region of physically allowed points {po,pi) 
is formed by a triangle po -l-pi < 1, Pi > 0. The blue convex 
region Q represents probability pairs achievable by mixtures 
of Gaussian states. The points lying in the light gray region 
indicate states that cannot be expressed as convex mixtures 
of Gaussian states. The green dashed line represents one of 
the non-Gaussianity witnesses W{a). (b) Detection scheme. 
Signal light beam impinges on a beam splitter (BS) and the 
outputs are detected by two APDs. Both single and coinci- 
dence rates are acquired by coincidence logic (CL). 



Pr] ^ G then also p ^ G- We can thus conservatively as- 
sume perfect detectors with unit efficiency and ii p-q ^ G 
is proven under this assumption, then it certainly holds 
also for p irrespective of the exact value of rj. 

In the experiment, the number of single detector clicks 
(RiA and Rib) as well as number of the coincidence chcks 
{R2) is measured for a given number of samples Rq of 
the state. Assuming perfect detectors with rj — I the 
vacuum-state fraction po is the probability that none of 
the detectors clicks, 



Po = 1 



RiA + Rib + R2 

Rn 



(3) 



The determination of pi is less trivial. We have 



n— 1 



i|^ = y:(i-T)>„, (4) 



n=l 



where T denotes the effective transmittance of the BS 
that also includes possible imbalance of the detection ef- 
ficiencies and other factors. Note that Ria and Rib de- 
pend on the whole photon number distribution p„, not 
just on Po and pi. We can nevertheless construct the 
following estimator. 



Pi, est 



RiA + Rib T2 + (1 - Tf R2 



Ro 2T{1 -T) Ra 

With the help of Eqs. ([3]) and ^ one can show that 



(5) 



Pl,CSt = Pi 



00 

E 

n— 3 



Pn- 



T2 - T" -h (1 - Tf - (1 - TY' 



2T{1 - T) 



hence pi,ost < Pi- Note that the term proportional to p2 
is absent in pi,est, so for rapidly decaying distributions 
the error is of the order of pa . With this lower bound on 
Pi at hand, the above criterion is still applicable, because 
apo -l-pi,cst > Wcia) implies that apo -I- pi > Wcia) as 
well. 

The estimation of pi is influenced by the effective im- 
balance of the detection channels T : (1 — T). Without 
loss of generality we can assume that T > i. It follows 
from Eq. ([5]) that pi^est decreases with increasing ratio 
T : (1 — T). Hence we should avoid underestimation of 
T which would result in overestimation of pi . An upper 
bound on T is provided by the ratio of single detector 
chcks, 



T,, 



RiA 



R 



lA 



Rib 



(6) 



We can thus 



It can be shown that T < Test for T > ^. 
safely use Test as a conservative estimate of T. 

Multimode witness: Many single-photon sources do not 
emit photons strictly into a single spatial and temporal 
mode. Let us therefore briefly sketch a proof [30| of the 
applicability of the witness to a generic multi-mode case. 



Let N denote the total number of modes involved and 
we define the total photon number n = X^j^i '^i' where 
rij is the number of photons in j-th mode. The esti- 
mated probabilities p„ then correspond to probability of 
no photon (n = 0) or one photon in total (n = 1) in 
the signal beam. Even in this multimode case, the max- 
imum of W(a) = apo + pi over all mixtures of Gaussian 
states is attained by a pure A^-mode Gaussian state. Any 
iV-mode pure Gaussian state can be prepared by combin- 
ing N single-mode squeezed states in a network of beam 
splitters. Moreover, the passive linear network described 
by a unitary matrix Un does not change the statistics of 
total photon number because UjynUj^ = n. It therefore 
sufHces to carry out the optimization over products of 
N pure single-mode Gaussian states which can be done 
analytically and one can prove that the bound on W{a) 
remains Waia) for arbitrary N. We can therefore ap- 
ply the witness also to multi-mode states without any 
limitation. 

Experimental setup: The experimental setup of a her- 
alded single photon source based on PDC is presented in 
Fig. [2j The detection part consists of three binary de- 
tectors (TR,SA,SB). The trigger detector (TR) yields a 
heralding output of the single photon PDC source. When 
a detection event is registered at this detector, we expect 
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FIG. 2: (Color online) Layout of the experimental setup. 
Continuous-wave frequency-multimode ultraviolet laser with 
central wavelength of 407 nm pumps a 2 mm long /3-barium 
borate (BBO) nonlinear crystal phase matched for type-2 de- 
generate PDC. Collinearly generated photon pairs with cen- 
tral wavelength of 814 nm are coUimated and separated from 
the pump beam by a dichroic mirror (DM) and spectrally 
limited by an interference filter (IF) to 10 nm. The orthog- 
onally polarized photons of a PDC pair are separated by a 
polarizing beam splitter (PBS) and coupled to single-mode 
optical fibers (SMF). Both outputs can be mixed with an at- 
tenuated (AT) infrared laser diode (LD) signal at fiber beam 
splitters (BS) to emulate dark counts of detectors and a noise 
component of the explored quantum state. The three out- 
put modes are then detected by binary single photon detec- 
tors based on sificon avalanche photodiodes (APD) operated 
in Geiger mode and actively quenched. The absolute quan- 
tum efficiency of the detectors is specified by a manufacturer 
to approximately 50%, while their relative efiiciencies were 
precisely measured prior to the experiment and found to be 
100%, 91.7 ± 0.2%, 92.2 ± 0.2% for TR,SA,SB channels, re- 
spectively. Electronic dark counts of the detectors in coinci- 
dence basis were found to be completely negligible. 



TABLE I: Estimated probabilities po and pi, and the corre- 
sponding witness I\W are shown for several different pump 
powers P and IF widths w (— denotes no filter). 



P[mW] 


w [nm] 


PO 


Pi 


AWixlO"*^] 


50 


2 


0.9124 


0.0875 


412 ±1 


50 


10 


0.8589 


0.1410 


1666 ± 3 


20 


10 


0.8425 


0.1574 


2370 ± 2 


50 


- 


0.7095 


0.2901 


14252 ± 17 


5 


- 


0.7296 


0.2704 


11825 ±15 



that an approximate single photon state is prepared in 
signal mode. The signal is divided by the beam split- 
ter (BS) with transmittance of Tbs = 0.522 ± 0.003 to 
the detection channels SA and SB. All single as well as 
two-fold and three-fold coincidence events between chan- 
nels TR, SA, and SB are registered by a fast coincidence 
logic unit. The overall splitting ratio T^^^ between the 
channels SA and SB has been conservatively estimated 
from the measured rates R\a and R\b using Eq. © and 
it agrees well with the independently measured Tbs £md 
relative detector efhciencies. 

The probabilities po a-nd p\ are estimated from the 
measured data using formulas ^ and ([5]). Due to condi- 
tioning on clicks of the trigger, i?o is given by the singles 
of the trigger detector, R\a and R^b by two-fold coinci- 
dences of TR&SA and TR&SB, respectively, while i?2 is 
actually given by the three-fold coincidences. The coinci- 
dence window is set to 2 ns. The results are summarized 
in Table I for different pump powers P and three differ- 
ent full width at half maximum (FWHM) of the interfer- 
ence filter IF (2 nm, 10 nm and without filter). Due to 
imperfect mode matching, in-coupling losses and ineffi- 
cient detectors, the vacuum term po dominates while the 
single-photon fraction is below 30% for all data shown. 
The contribution of higher photon terms is very small, 
\ — Pq— P\ < lO"'*, so the generated state is very close 
to an attenuated single photon. In the experiment, this 
is indicated by a very small ratio of three-fold to two- 
fold coincidence rates (less than 10"'^). For example, 
for P = 50 mW and ui = 10 nm we have i?2 = 605, 
RiA = 1-259 X 10^ and Rib = 1-192 x 10^ per 100 s. The 
statistical uncertainty of the estimated po and p\ deter- 
mined assuming Poissonian statistics is less than 2 x 10~ 
for all data shown (one standard deviation) . Due to very 
low three- fold coincidences, the sumpo+Pi exhibits much 
lower statistical uncertainty than the difference po — pi- 

Results: We have verified that the generated states 
cannot be expressed as a mixture of Gaussian states by 
using the non-Gaussianity witness. For each data set 
{po,Pi) we have calculated AW = apo+pi — Wo{a) and 
maximized the difference over all a. The resulting maxi- 
mal AW are listed in Table I. We can see that AW > 
in all cases and the bound Wcia.) is always surpassed 



TABLE II: The same as Table I, nrd indicates the amount 
of noise injected from LDl and LD2 into trigger and signal 
detectors, P — 50 mW, and IF width w — 10 nm. 



Wrcl 


Po 


Pi 


f^opt 


AW^[xl0-«] 


0.0 


0.8195 


0.1804 


0.94018 


3479 ± 7 


0.1 


0.9073 


0.0926 


0.98389 


406 ±3 


0.2 


0.9408 


0.0591 


0.99332 


42 ±2 


1.0 


0.9777 


0.0222 


0.99903 


-84 ±1 



by many standard deviations. Next, we investigate the 
influence of background noise on the source properties. 
For this purpose we inject light from laser diodes LDl 
and LD2 into trigger and signal detection blocks, respec- 
tively. Noise from LD2 emulates background noise of the 
source while noise coming from LDl effectively increases 
dark count rate of the trigger thus increasing the vacuum 
fraction po- Table II shows the results obtained when the 
amount of injected noise is the same in both blocks and 
rirei indicates the normalized noise strength. With in- 
creasing noise we can clearly observe transition to the 
regime where AW < 0, as also illustrated in Fig. 3. 

Discussion: Let us briefly compare our results with 
other non-classicality measures. Since po > 0.5 for all 
the measured states, their Wigner function is always pos- 
itive in the origin, 1^(0) = ^((-1)") > ^^ > 0, 
where it is expected to exhibit maximum negativity 
W{0) = —- for a pure single photon state [3l[. On 
the other hand, all the measured states cannot be ex- 
pressed as a mixture of coherent states, therefore they 
are non-classical [32|. The non-classicality can be quan- 
tified by a g(^) parameter defined as g^'^\0) = "t^t;^ 
[33| . Sub-Poissonian photon number statistics is in- 



dicated by g^'^\0) < 1. The state produced by our 
source can be excellently approximated by a density ma- 
trix pT = Po\0){0\ +pi|l)(l| + (1 -Po -Pi)|2)(2| be- 
cause higher photon terms are exponentially suppressed 
due to very low parametric gain in the nonlinear crys- 



tal. In this case we find g^^''(0) 



2(l-pi-po) 
l2(l-po)-pil^ 



yielding 

^'^•'(O) < 0.3661 for all states. Simultaneously, all the 
results exhibit very strong photon anti-correlation effect, 
witnessed hy a — „^°?,^ < 0.3706 which violates the 
classical inequality a > 1 [15|. The limits are given by 
data for rij-d = 1 in Table II. All the above parameters are 
monotonously decreasing as less noise is imposed by LDl 
and LD2 and for nrci = 0.2 we already have a = 0.0521 
and g(2)(o) = 0.0519. 

In conclusion, we have examined a source producing 
approximate single-photon states with positive Wigner 
function but exhibiting strong photon anti-correlation 
and we have unambiguously proved that the generated 
states cannot be expressed as mixtures of Gaussian 
states. In comparison to the witness based on negativ- 
ity of the Wigner function [3l| , the present criterion can 
identify a high nonclassicality of a much wider class of 
single photon sources, while avoiding the need for de- 
manding estimation of complete photon number distribu- 
tion or complicated data processing |3J]. Consequently, 
the presented criterion is particularly useful for evalua- 
tion of single-photon sources where negativity of Wigner 
function cannot be observed [35|. 
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FIG. 3: (Color online) The multiphoton contribution 1 — po — 
pi is plotted as a function of pi . Symbols indicate experimen- 
tal data for several different levels of noise added simultane- 
ously by LDl and LD2 (red squares) and by LDl only (green 
circles). All other parameters were fixed, P = 50 mW and 
w = 10 nm. Error bars stand for three standard deviations, 
horizontal error bars of pi are smaller than the size of the 
symbols. The solid blue curve represents the boundary given 
by Eq. ([l} , all points on the right of this curve correspond to 
states that cannot be obtained as mixtures of Gaussian states. 
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